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^ : ABSTRACT 

C ' Context. Transition disks typically appear in resolved millimeter observations as giant dust rings surrounding their young host stars. 

^~r^' More accurate observations with ALMA have shown several of these rings to be in fact asymmetric: they have lopsided shapes. It 

■^^ I has been speculated that these rings act as dust traps, which would make them important laboratories for studying planet foiTnation. 

It has been shown that an elongated giant vortex produced in a disk with a strong viscosity jump strikingly resembles the observed 

asymmetric rings. 

Aims. We aim to study a similar behavior for a disk in which a giant planet is embedded. However, a giant planet can induce two 
Q, ] kinds of asymmetries: (1) a giant vortex, and (2) an eccentric disk. We studied under which conditions each of these can appear, and 

how one can observationally distinguish between them. This is important because only a vortex can trap particles both radially and 

azimuthally, while the eccentric ring can only trap particles in radial direction. 
(~ l ' Methods. We used the FARGO code to conduct the hydro-simulations. We set up a disk with an embedded giant planet and took a 

f^ I radial grid spanning from 0.1 to 7 times the planet semi-major axis. We ran the simulations with various viscosity values and planet 

masses for 1000 planet orbits to allow a fully developed vortex or disk eccentricity. Afterwards, we compared the dust distribution in 

a vortex-holding disk with an eccentric disk using dust simulations. 

Results. We find that vorticity and eccentricity are distinguishable by looking at the azimuthal contrast of the dust density. While 
^ I vortices, as particle traps, produce very pronounced azimuthal asymmetries, eccentric features are not able to accumulate millimeter 

dust particles in azimuthal direction, and therefore the asymmetries are expected to be modest. 
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fT^ ■ 1. Introduction Bimstiel et al. 2013). For instance, Regaly et al. (2012) and 

t^^ ' Lyra & Mac Low (2012) have shown that if gas accumulates at 

■ Most models of protoplanetary disks that aimed at fitting obser- ^^^^ ^^^^^-^^ -^ ^^^ ^^^^^ ^^^j^ ^^ ^ sharp jump in viscous stress 

l''^T'J'^!'''''^'^'^T^ t^- ^.^^f,° u ^^^' ^'"^'^ ^^ , or resistivity, the gas "bump" can become Rossby-unstable and 
2006; Woitke et al. 2009; Mm et al. 201 1) because until recendy, ^^^^^^ ^J anti-cyclonic vortices. Initially, a number of small 
^ , with a few exceptions, observations did not have enough spa- ^^^-^^^ ^^^ generated in these models that eventually merge to 
tial resolution and sensitivity to detect strong deviations from ^^^^ ^ -^^^ crescent-shaped vortex, or in other words, a lop- 
axial symmetry in these disks. Recent observations, however, ^-^^^ ^. ^^^^^ ^^^.^-^^^ ^^^j^ ^j^^ ^^^.^ ^^ ^^^ ^^ ^^ ^ 
have unveiled some non-axisymmetric structures in protoplane- educed by a massive planet. As the planet opens up a gap in 
tary disks and have changed the situation. Scattered-light images ^^^ ^-^^ ^^^ ^^^^ ^^^ ^^ j^^^^ hampers) the gas flowing to- 
_ of a number of disks around Herbig Ae/Be stars have revealed ^^^^ ^j^^ ^^^^^ ^j^^ accumulates at the outer edge of the gap, 
^ • S^^f ^^ structures such as spirals and rings (e.g^Fukagawa et al. ^^^^^^, ^n^table, and produces vortices (Koller et al. 2003; Li 
■ ■ 2004; Oppenheimer et al. 2008; Muto et al. 2012). Millimeter ^^ ^j 2005; De Val-Borro et al. 2007; Lin & Papaloizou 2011; 
continuum maps show asymmetries as well (Pietu et al. 2005). Lin 2012). However, it is less clear how long these vortices can 
Spatially resolved mm maps consistently show that transitional c.j^ive ("Meheut et al 2012") 
disks -protoplanetary disks with a strong deficiency of dust in- 
ward of a few AU, or in other words, a disk with a huge in- 
ner hole (Andrews et al. 201 1)- appear as rings on the sky of- In addition, a massive enough planet can produce an eccen- 
ten with a lumpy structure along the ring (Brown et al. 2009; trie gap and thus generate another kind of large asymmetry. Kley 
Casassus et al. 2013). Near-infrared scattered-light images of & Dirksen (2006) showed that a planet with a mass larger than 
some of these disks also clearly show non-axisymmetric struc- 3Mjup in a disk with a viscosity ff ^ 4 x 10""' is able to produce 
tures (Rameau et al. 2012). Altogether, non-axisymmetry turns an eccentric disk which means that the gas parcels in the disk are 
out to be a regular phenomenon in disks, particularly in transi- on elliptic orbits. Because of Kepler's laws, the gas velocities at 
tion disks. the apocenter of the elliptic orbits are lower than at the pericen- 
Recent theoretical work has aimed at explaining the ori- ter As a result of this "traffic jam", the density p is expected to 
gin and physical processes involved in these asymmetries (e.g. be higher at apocenter than at pericenter for a nearly constant 
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mass flux {pu^) along the ellipse. This also leads to a lopsided 
disk structure. 

The question arising here is how the basic diff'erences be- 
tween two models, vortex and eccentric disk, can aff'ect the ob- 
servational features of transition disks. One important difference 
between the eccentric structure and vortex is the density contrast. 
A higher contrast is achievable if the planet produces a more ec- 
centric disk, in which case the eccentricity of the orbit presum- 
ably becomes apparent in the image. Another notable difference 
is that the eccentric traffic jam is not comoving with the fluid, but 
instead is nearly stationary in the inertial frame (Kley & Dirksen 
2006). In simple words, similar to galactic spiral arms, the mass 
flows in from one side and flows out at the other side. As a re- 
sult, this kind of asymmetry cannot trap particles, and the dust 
asymmetry would be similar to the gas. This has been confirmed 
by Hsieh & Gu (2012), who showed that in an eccentric disk 
with a ~ 5Mjup planet, the dust density follows the gas density 
pattern. Instead, anticyclonic vortices, which are local pressure 
"bumps", act as dust traps (e.g Barge & Sommeria 1995; Klahr 
& Henning 1997; Lyra et al. 2009). The dust trapping may lead 
to an even stronger azimuthal asymmetry in the dust continuum 
(Birnstieletal. 2013). 

The objective of this paper is to study the observational ap- 
pearance of asymmetries caused by a planet and answer the 
question when we observe a transition disk and see a large-scale 
asymmetric ring, whether this is due to a vortex or an eccentric 
disk. We aim to find out under which conditions eccentricities or 
vortices are the dominant sources of azimuthal asymmetries. 



2. Method and simulations 

Both a vortex and an eccentric disk, which cause an asymmetry 
in a gas disk, can be produced by a massive planet and are af- 
fected by the physical parameters of the disk in the same manner, 
although they have diff'erent formation mechanisms. In Sec. 2.1 
we briefly diskuss the formation of these asymmetries and the 
parameters that can aff'ect them, then we describe our numerical 
methods. 



2.1. Vortex and eccentric shape 

A large-scale vortex can be created by a Rossby-wave instabil- 
ity (RWI) and live long under suitable conditions. A RWI can 
be considered as the rotational analogy of the Kelvin-Helmholtz 
instability (Lovelace et al. 1999). Li et al. (2000) studied the 
condition of RWI formation in a disk with either a density bump 
or a density jump and concluded that a 10-20% surface density 
change over a radial length scale of about the thickness of the 
disk is enough to perturb the rotational velocity of the gas and 
produce the RWI. Wherever this condition is satisfied in the disk, 
the instability occurs. De Val-Borro et al. (2007) showed that 
not only planets as massive as Neptune and Jupiter are able to 
excite the RWI and generate vortices on both edges of the gap, 
but the vortices can also live longer than 100 orbits. A 3D study 
by Meheut et al. (2012) showed that RWI can survive for long 
timescales (~ hundred orbits), and the authors explained that, 
as long as the overdensity is sustained permanently by an exter- 
nal driving force such as (in our case) a planet, the vortex does 
not decay and is not destroyed by mechanisms such as elliptical 
instability. Therefore, parameters such as viscosity and planet 
mass, which control the accretion rate through the gap and con- 
sequently the density gradient in the outer edge of the gap, can 
influence the formation and the lifetime of the vortex. 



An eccentric gap is another by-product of gap formation by 
a massive planet and is affected by parameters similar to those 
that influence a vortex. Kley & Dirksen (2006) investigated the 
response of disk eccentricity by changing planet mass and phys- 
ical parameters of the disk such as viscosity and temperature. 
They showed that if the planet is massive enough to clear a gap 
that reaches the outer 1:2 Lindblad resonance, the disk becomes 
eccentric. Like vortex formation, the eccentric structure is also 
altered by viscosity and planet mass. High viscosity and lower 
planet mass narrow the gap and reduce the eccentricity of the 
disk. We ran 132 hydrodynamical simulations to study how these 
parameters affect the formation and lifetime of a vortex and/or 
an eccentric feature. 



2.2. Hydro-simulations 

We used the locally isothermal version of the FARGO code 
(Masset 2000). Our basic model is a 2D viscous disk with a mas- 
sive planet that opens a gap. We altered viscosity and mass of the 
planet to see how they affect the asymmetric features. We per- 
formed 12 runs in inertial frame for 1000 orbits. To distinguish 
between a vortex and an eccentric shape, we followed the disk 
surface density evolution at a high time-resolution. To do this, 
we conducted an extra set of ten runs per model to extend every 
output in each simulation for another two orbits, which saved 
100 outputs during each orbit. 

We considered a flared disk with an aspect ratio of h = 
Q.05{r/rpf'-^^. Considering the aspect ratio definition h{r) = 
H/r - cJvk, where c, is the local sound speed and \'k is the 
Keplerian velocity, the temperature profile scales as T oc r '^^. 
The surface density follows the relation E = So {rlvpY^, with 
So = 2 X 10"^ Mi,lrl,, and to keep the disk viscously stable, we 
used a viscosity of v = aCsH (Shakura & Sunyaev 1973). The 
viscosity parameter a had the values 10"^, lO"-', and 10"^ in our 
models. The disk extended from rn,in = O.lr,, to r^^^ - l.Orp, 
where r^ is the orbital radius of the planet and was used as 
length scale. The disk is covered by A^^ x Ns = 512 x 757 grid 
cells with logarithmic radial spacing. We chose this resolution 
to avoid numerical problems caused by the high density gradient 
at the gap edge and to have squared cells. We used the FARGO 
non-reflecting boundary condition to reduce the effect of wave 
reflections from the boundaries. 

The planet mass in our simulations had the values Mp = 
5, 10, 15, and 20 Mjup with Mjup = lO^^M*. The planet was held 
at r^ = 1 in a circular orbit. The potential of the planet <p was 
softened by the parameter e = Q.6Rh to avoid a singularity; 



GMp 

'(r^ + e^y 



(1) 



where Rh is the planet Hill radius of the planet and G is the 
gravitational constant. 

To apply our results to the observed features of the transition 
disks, we need the planet orbital radius and mass of the star in 
physical units. The remaining quantities can be calculated using 
these two numbers. In this work we considered a planet that or- 
bits a solar-mass star at a distance of rp = 20AU. Therefore, the 
disk extends from r,nin = 2AU to r^ax = 140AU and has a mass 
of Mdisk = O.OO8M0. The initial surface density So is therefore 
So ~ 4.44 (r/20AU)-'g cm^^. 
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Fig. 1. Evolution of a vortex-holding disk and an eccentric disk 
from 1000 orbits (left column) to 1002 orbits (right column). The 
top panels show how a vortex moves during two orbits of the 
planet. The bottom panels show the same snapshots as the top 
panels, but for an eccentric disk. Clearly, the eccentric feature 
is fixed during the two planetary orbits. The crosses display the 
planet position. 



2.3. Recognizing a vortex from an eccentric shape 

We distinguished a vortex from an eccentric asymmetry in our 
hydro-simulations by checking whether the structure was fixed 
or co-moving with the gas in the inertial frame. While a vor- 
tex revolves around the star with the local orbital frequency, an 
eccentric feature is almost fixed with small precession (Kley & 
Dirksen 2006). Accordingly, we followed the structures in all 
models to decide whether the asymmetry was a vortex or an ec- 
centric feature (Fig. 1). 



2.4. Dust simulations 

To compare the dust distribution between a vortex and an ec- 
centric feature, we conducted dust simulations for the models in 
Fig. 1. We used the numerical code developed by Zsom et al. 
(201 1), which solves the restricted three-body problem for dust 
particles in 2D spherical coordinates, taking into account gas 
drag that results in radial drift (for details see Paardekooper 
2007; Zsom et al. 2011). To mimic the turbulence mixing, we 
allowed the particles to pace every time step (At) in a random 
direction with a step size of / - -sjDpAt, where Dp - v/(l + t^^) 
is the diffusion coefficient for the particles with dimensionless 
stopping time Tj (see Youdin & Lithwick 2007). The viscosity 
V was taken to be self-consistent with our hydrodynamical sim- 
ulations. Because we are only interested in the dust dynamics, 
we did not consider the dust coagulation and fragmentation. We 
ran two sets of simulations for two particle masses of 0.05 and 
0.5g. These masses correspond to sizes ~ 2 and ~ 5mm in our 
simulations which suffer the highest and moderate radial drift 
at ~ 5QAU. We distributed 10000 dust particles with the same 
density profile as the gas at 1000 orbits and allowed particles to 
evolve for 200 orbits. 
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Fig. 2. Dependency of eccentricity and vortex lifetimes on planet 
mass and viscosity. The color represents the lifetime of vortices, 
and the eccentric disks are marked by a dot. For Mp - 15Mjup 
and a - 10"^, the asymmetries are extremely mixed and indis- 
tinguishable after 600 orbits. 



3. Results 

Figure 2 summarizes of the results of our models in terms of 
the presence or absence of a vortex and/or an eccentricity. The 
horizontal and vertical axes are the mass of the planet and the 
logarithm of the a viscous parameter According to Fig. 2, an 
eccentric shape is a very common feature in our models. This 
figure shows that a vortex can be created and survive until the 
end of the simulations for a low viscosity, i.e., a = lO"''. 

In none of the simulations with the highest viscosity param- 
eter a = 10"^, was a vortex formed. In the model with a 5Mjup 
planet, a narrow circular gap was created and stayed until the end 
of the simulations. A more massive planet Mp - lOMjup starts 
to form an eccentric gap from 200 orbits and establishes this af- 
ter 400 orbits (~ 3.6 x 10'* years). For more massive planets, the 
disk becomes eccentric even at 200 orbits (~ 1 .8 x lO'* years). 

For an intermediate value of the viscosity, i.e., a = 10"^, the 
disk passed through three stages. First a vortex appeared in the 
disk and survived until 200-300 orbits. The lifetime of the vor- 
tex greatly depends on mass of the planet. The more massive the 
planet, the longer the vortex survives. During the second stage, 
which lasts until orbit 800 for Mp - 5Mjup and until orbit 500 for 
the remaining masses, the features are mostly indistinguishable 
and mixed. Therefore, even if we observe a moving denser re- 
gion or a vague asymmetry, we consider it to be neither a vortex 
nor an eccentric disk. Eventually, the planet succeeds in reform- 
ing the gap into an eccentric shape and the disk keeps this until 
the end of the simulation. 

The models with a low eccentricity of a = lO"'* produce 
the most interesting results for the vortex formation. For lower 
planet masses, i.e., Mp - 5Mjup and lOMjup, the vortex is the 
dominant feature. Although the disk is slightly eccentric, the 
asymmetric density feature caused by it is not strong enough to 
be distinguished, and therefore the vortex remains the strongest 
asymmetry in the disk. On the other hand, the eccentric feature in 
the models with more massive planets Mp = 15Mjup and 20Mjup 
is strong to such a degree that it can mix with the vortex. When 
the two features overlap, the density contrast is the the highest. 

Figure 3 illustrates the results of the dust simulations for the 
two disk models in Fig. 1 . The first row of Fig. 3 shows the dust 
distribution for the disk with the vortex. In this case, particles 
move toward the vortex, and after 200 orbits most of the dust is 
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Fig. 3. Dust distribution comparison at 1019, 1051, 1102, and 
1197 orbits. We chose these snapshots because the vortex is 
located at the same azimuthal position of the eccentric feature 
which facilitates the comparison of the azimuthal particle con- 
centration. The top and middle rows reflect the diff'erence be- 
tween the dust trapping in a vortex-holding disk and an eccentric 
disk. The light and dark colors represent the 0.05g and 0.5g par- 
ticles. In the bottom panels, we plot the dust enhancement factor 
against azimuth at each corresponding snapshot. 



accumulated in the vortex center. The middle panels of Fig. 3 
display the the dust distribution in the eccentric case. In contrast 
to the vortex case, the particles are only trapped in the radial 
direction, which produces a ring that is wide because of the high 
diffusion (a = 10"^). In the top panels the turbulence stirring is 
not efficient because of the low viscosity and the results would be 
very similar to those of the models without stirring. However, the 
particles in middle panels are less concentrated around the outer 
edge of the gap because of the higher diffusion coeflicient. After 
the semi-steady state is reached, no strong azimuthal asymmetry 
is produced. We define the dust enhancement factor f as 



/W = 



(Mdust/Mgas), 



(A/dust/A/gas)f=1000 orbits 



(2) 



would happen for both cases (see e.g Pinilla et al. 2012), while 
particle trapping in azimuthal direction can only happen in the 
vortex case. Because a traffic jam in an eccentric disk is no par- 
ticle trap, the azimuthal contrast in the continuum would be as 
high as the contrast in the gas surface density, which can reach 
values of S/Sq ~ 3 (Fig. 1), and the dust azimuthal extension 
of the asymmetry would be similar to gas. Conversely, vortices 
are indeed particle traps and can be long-lived (up to 10^ years) 
when very massive planets (15-20 Mjup) and moderate viscosity 
(a = 10 "^) are able to maintain the density bumps and the RWl. 
Thus, the vortices can create strong variations in the azimuthal 
dust distribution of millimeter particles and the asymmetry fea- 
ture would be less azimuthally extended for the dust than for the 
gas (Bimstiel et al. 2013). Hence, stronger and more compact 
dust azimuthal asymmetry is expected for a vortex than for an 
eccentric disk. The results shown in Fig. 3 give us a good im- 
pression of what is expected with millimeter observations with 
high sensitivity and angular resolution. It is important to notice 
that neither in a vortex nor in an eccentric case, we expect the 
location of the planet along its orbit to be correlated with the 
location of the vortex or elhpse of disks. 

5. Conclusion 

One exciting explanation of the wide gaps observed in transi- 
tion disks is the interaction of a massive planet with the disk. 
When the planet is massive enough and moderate values for the 
disk viscosity are assumed, two different scenarios can create 
azimuthal variations in the gas surface density: eccentric disks 
and a vortex excited by the RWL The azimuthal contrast of the 
gas surface density varies for both cases and is much higher for 
the vortex situation. In this case, particle trapping is possible and 
very high azimuthal variations are expected at millimeter wave- 
lengths, as is currently observed with ALMA. 
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which is calculated for r e [30, 140]AU and the azimuthal slice 
between (p and (f> + Acp. We plot f{(p) in the bottom panels of 
Fig. 3 to compare the dust azimuthal accumulation of our mod- 
els. These plots show that the dust enhancement factor in a vor- 
tex is much higher than the value for the eccentric case. This 
means that, the vortex produces a stronger azimuthal asymmetry 
in a narrower azimuthal range than the eccentric disk. 

4. Observational consequences 

The difference between azimuthally density contrast in a vortex- 
holding and an eccentric disk is a very useful tool for rec- 
ognizing them observationally. Our results show that we need 
~ 50 X (/? J /20A {/)'■'' years {Rs is the asymmetry position) to 
distinguish an eccentric feature from a vortex using the compar- 
ison of the position of the azimuthal density variation. However, 
different observational features at (sub-)millimeter wavelengths 
are fortunately diskernible in one single snapshot. When a planet 
carves out a gap in a disk, particle trapping in the radial direction 
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